************************* LW^st siiare Lattice ***************************** 
************************ Noise cancelling ****************************/ 

/* Example for constant saturation approach to noise cancelling */ 
/define LAMBDA 0.95 


int 

resets 

int 

passes. 

int 

Bat_f actor. 

int 

*6ignal_l. 

int 

*6ignal"2. 

int 

♦target 1, 

int 

♦target 2) { 


int i, ii, k, in, n, contractions- 
static int *6 a, *6_b, *out_a, *out_b; 
static float DeTta_sqr, scale, noise_ref; 

if ( reset « TRUE) { 
s_a « Bignal_l; 
s_b = signal_2; 
out_a « target_l; 
out"b «= target_2; 
scale «= 1.0 /4160.0; ' * 

/* noise canceller initialization at time t«0 */ 

nc[0] .berr = 0*0; 
nc[0] .Gamma =1.0; 

for(m«0; m<NC_CELLS; m++) { 
nc(in].err_a = O.O; 
nc[in].err^b «= 0.0; 
nc(m].Roh_a « 0.0; 
nc(m).Roh_b = 0.0; 
nc(in]. Delta = 0.0; 
ncjoj .Fswsqr = 0.00001; 
nc[m].Bswsqr = 0.00001; 

} 

} 

for(k«=0; k<passes; k++) { 

contraction » FALSE; ^ ^ ^ ^ , . 

for(m«0; m< NC^CELLS; m++) { /* Update delay elements */ 

nc[m].berrl «= nc[m)-berr; 

nc(m) .Bswsqrl - nc [m] . Bswsqr; 

} 

noise ref « sat_f actor * log (1.0 - (*s_a) * scale) 

- log(l*0 - (*s_b) * scale) ; 
nc(0],err_a = log(l,0 - (*s_a) * scale); 
nc(0].err_b = log(1.0 - (*s_b) * scale); 


++S a; 


• I 


++s_b ; 

nc[0] . f err 
nc[ 0] . berr 
nc[0] • Fswsqr 
nc(0 j .Bswsqr 


noise_ref ; 
noise_ref ; 

LAMBDA * nc[0]. Fswsqr + noise_ref * noise_ref; 
nc(0] . Fswsqr; 


/* Order Update */ 

for(n-=l;( n < NC_CELLS) (contraction 

/* Adaptive Lattice Section */ 


FALSE) ; n-M-) { 


ii- 


n-l; 
n-l; 


nc(m). Delta *« 
nc(m] .Delta +«= 
Delta_Bqr = 


nc[n] . fref 
nc(n] .bref 

nc(n] . f err 
nc(n j .berr 


LAMBDA; 

nc(mj.bem * nc(iii].ferr / nc[ni] •Gamma 
nc[m]«Delta * nc(m]. Delta; 

-nc[m]. Delta / nc(m] .Bswsqr 1; 
-nc ( m ] . De It a / nc ( m ] .Fswsqr ; 


= nc(m].ferr + nc[n].fref * nc[m].berrl; 
= nc(in].berrl + nc(n].bref * nc[m].ferr; 


nc(n]. Fswsqr = nc(m]. Fswsqr - Delta^sqr / nctn] .Bswsqrl, 
nc(n). Bswsqr = nc(ia).Bswsqrl - Delta_sqr / nc[m] . Fswsqr ; 


if{ (nc[n] .Fswsqr + nc[h) .Bswsqr) > O.OOOOl 
nc(n]. Gamma = nctm]. Gamma - nc(m].berrl * 
if (nc[n] .Gamma < 0.05) nc[n]. Gamma = 0.05; 
if (ncfnj .Gamma > 1.00) nc(n]. Gamma =1.00; 

/* Joint Process Estimation Section */ 


II (n < 5) ) { 

nc[m].berrl / nc[m] . Bswsqrl ; 


nc(m].Roh_a *= 

nc(m] .Roh_a += 

nc(m].)c_a = 

nc(n].err_a = 

nc[m].Roh_b *= 

nc[m] .Roh_b += 

nc[mj.k_b = 

nctn].err_b = 


LAMBDA; 

nc[m].berr * nc(m].err_a / nc[ra]. Gamma 
nc(mi.Roh_a / nc(m] . Bswsqr ; 
nc[m].err_a - nc(m].}c_a * nc(m]-berr; 


LAMBDA; 

nc(m].berr * nc(m].err_b 


nc(m] .Roh_b 
nc(m] .err_b 


. . _ / nc[m] .Gamma 
/ nc[m] .Bswsqr; 
- nc(ml.k_b * nc[m).berr; 


} 

else { 

contraction «= TRUE; 

for ( i=n ; i<NC_CELLS ; i++) { 


nc[i] .err_a 
nc [ i ] . Roh_a 
nc[i] .err_b 
nc( i] ,Roh_b 
nc[i] .Delta 
nc[i) .Fswsqr 
nc[ i] . Bswsqr 
nc[ i] . Bswsqrl 


0.0; 
0,0; 
0.0; 
0.0; 
0.0; 

0. 00001; 
0.00001; 
0.00001; 


} 

} 

) 

*out_a++ » (int) ( (-exp(nc(ii] .err_a) +1.0) / scale) ; 

*out_c++ o (int) ( (-exp(nc[ii] .err_b) +1.0) / scale) ; 



Least Square Lattice *************************** 


